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Context. Plaskett's star (HD 47 129) is a very massive O + O binary that belongs to the Mon OB 2 association. Previous work suggested 
that this system displays the Struve-Sahade effect although the measurements of the secondary radial velocities are very difficult and 
give controversial results. Both components have powerful stellar winds that collide, and produce a strong X-ray emission. 
Aims. Our aim is to study in detail the physical parameters of this system, and to investigate the relation between its spectral properties 
and its evolutionary status. 

Methods. We present here the analysis of an extensive set of high resolution optical spectra of HD 47 129. We use a disentangling 
method to separate the individual spectra of each star. We derive a new orbital solution and discuss the spectral classification of both 
components. A Doppler tomography technique applied to the emission lines Ha and He n A 4686 yields a Doppler map that illustrates 
the wind interactions in the system. Finally, an atmosphere code is used to determine the different chemical abundances of the system 
components and the wind parameters. 

Results. HD 47 129 appears to be an 08 III/I + 07.5 III binary system in a post RLOF evolutionary stage, where matter has been 
transferred from the primary to the secondary star. The He overabundance of the secondary supports this scenario. In addition, 
the N overabundance and C underabundance of the primary component confirm previous results based on X-ray spectroscopy and 
indicate that the primary is an evolved massive star. We also determine a new orbital solution, with Mp sin 3 i = 45.4 ± 2.4 M Q and 
M s sin 3 i = 47.3 ± 0.3 M G . Furthermore, the secondary star has a large rotational velocity (v sin i ~ 300 km s~') that deforms its 
surface, leading to a non-uniform distribution in effective temperature. This could explain the variations in the equivalent widths of 
the secondary lines with phase. We suggest that the wind of the secondary star is confined near the equatorial plane because of its 
high rotational velocity, affecting the ram pressure equilibrium in the wind interaction zone. 

Key words, stars: individual: HD47129 - binaries: spectroscopic - stars: fundamental parameters 



1. Introduction 

Plaskett's star (HD 47 129) is a non-eclipsing early type binary 
which was analysed for the first time by Plaskett (1922). He 
found doubled spectral lines and described the system as consist- 
ing of a bright 05e primary component and a fainter secondary 
component with weak and diffuse lines. When he determined 
the elements of the spectroscopic orbit, the large mass of this 
system attracted his attention ((M P + M s )sin 3 i - 138.9 M ). 
This result motivated Struve (1948) to re-observe this object in 
1947-1948 over two complete orbital cycles. He detected that 
some lines in the secondary spectrum (Hy, He i AA 4026, 4471) 
were weaker when this component was receding. This effect is 
nowadays called the Struve-Sahade effect (S-S effect hereafter) 
and is currently defined as the apparent strengthening of the sec- 
ondary spectrum of a binary system when the secondary is ap- 
proaching, and the corresponding weakening of its lines when 
it is receding. Struve (1948) also noted differences in the in- 
tensity of the secondary spectrum from one cycle to another, 
at a given orbital phase. He explained these variations with the 
presence of streams of gas in the system. A third observational 
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campaign, conducted between 1956 and 1958 by Struve et al. 
(1958), allowed these authors to conclude that the behavior of 
the radial velocities of the secondary component suggests that it 
ejects matter in an irregular manner, and that these ejections are 
larger over the hemisphere facing the primary star. Hutchings & 
Cowley (1976) confirmed the presence of the S-S effect in the 
system and suggested that the primary star has a stellar wind 
which is stronger on the rear side of the star with respect to its 
orbital motion. IUE data allowed Stickland (1987) to determine a 
radial velocity curve and an orbital solution for the primary com- 
ponent, but not to identify the secondary spectral lines. Bagnuolo 
et al. (1992) used a tomographic technique on the same IUE 
data to separate the two components of the system. They found a 
mass ratio ofl.l8±0.12 (the secondary being more massive than 
the primary) and spectral types of 07.5 1 and 06 1 for the primary 
and the secondary, respectively. Given that they estimate an in- 
clination around 71°, they found masses equal to 42.5 M for the 
primary star and 51.0 M for the secondary star. The secondary 
lines appeared rotationally broadened (v sin i - 3 10+20 km s~') 
but are visible throughout the orbit. Finally they found that the 
secondary contribution to the cross-correlation functions is rel- 
atively constant with orbital phase. Stickland (1997) also used 
IUE data and cross-correlation functions to study specifically 
the S-S effect in Plaskett's star. The conclusion was that the sec- 
ondary radial velocities "make no sense dynamically", and that 
they are probably due to the presence of gaseous streams or en- 
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Table 1. Journal of the observations of HD 47 129. The first column 
indicates the instrument used to obtain the data, the second lists the 
Heliocentric Julian Date (HJD), the third provides the orbital phase (O) 
as computed from our orbital solution (see Table 2) and the fourth and 
fifth ones give the radial velocities used to compute the orbital solution, 
respectively for the primary and secondary stars. The dispersion on the 
radial velocities are on average less than 2.0km s" 1 for the primary, and 
of the order of 20 km s _1 for the secondary component. 

Instrument HJD O RV(kms _1 ) 

- 2 450 000 Primary Secondary 



FEROS 


2335.565 


0.28 


230.8 


-149.2 


FEROS 


2336.532 


0.35 


197.7 


-141.0 


FEROS 


2336.551 


0.35 


192.9 


-119.7 


FEROS 


2337.527 


0.42 


120.6 


-117.8 


FEROS 


2337.534 


0.42 


118.3 


-118.2 


FEROS 


2338.528 


0.49 


60.8 


-29.2 


FEROS 


2339.535 


0.56 


-28.4 


- 


AURELIE 


3648.675 


0.49 


- 


- 


AURELIE 


3650.662 


0.63 


- 


- 


AURELIE 


3652.673 


0.77 


- 


- 


AURELIE 


3654.663 


0.91 


- 


- 


FEROS 


3738.692 


0.74 


-169.0 


134.5 


FEROS 


3772.629 


0.10 


153.9 


-151.4 


FEROS 


3772.640 


0.10 


158.5 


-129.7 


AURELIE 


3775.486 


0.30 


- 


- 


AURELIE 


3776.486 


0.37 


- 


- 


AURELIE 


3777.478 


0.44 


- 


- 


AURELIE 


3778.462 


0.51 


- 


- 


AURELIE 


3779.318 


0.57 


- 


- 


AURELIE 


3780.288 


0.63 


- 


- 


FEROS 


3796.515 


0.76 


-171.7 


124.9 


FEROS 


3797.521 


0.83 


-152.0 


145.2 


FEROS 


3798.510 


0.90 


-80.9 


147.1 


FEROS 


3799.511 


0.97 


-23.2 


74.6 


FEROS 


3800.507 


0.04 


74.5 


- 


AURELIE 


3843.329 


0.01 


- 


- 


AURELIE 


4033.700 


0.24 


- 


- 


AURELIE 


4034.692 


0.30 


- 


- 


AURELIE 


4035.672 


0.37 


- 


- 


SOPHIE 


4162.432 


0.18 


205.9 


-181.6 


SOPHIE 


4164.307 


0.31 


219.7 


-190.5 


SOPHIE 


4168.468 


0.60 


-78.6 


155.3 


SOPHIE 


4172.305 


0.86 


-123.5 


187.6 


SOPHIE 


4174.333 


0.00 


31.1 


- 


SOPHIE 


4176.327 


0.14 


208.5 


-135.2 


SOPHIE 


4182.321 


0.56 


-45.3 


100.2 


SOPHIE 


4184.350 


0.70 


-164.5 


160.0 



velopes in the system. However, Bagnuolo et al. (1999) used the 
same method and concluded that there is no evidence of the S-S 
effect in HD 47 129. 

Both components of Plaskett's star are very massive and 
luminous O-stars, and they thus both possess strong stellar 
winds that collide, producing X-ray emission in the interaction 
zone. An X-ray study of this object has been made by Linder 
et al. (2006), using new XMM-Newton data and ROSAT archival 
data. The RGS (Reflection Grating Spectrometers of the XMM- 
Newton satellite) data showed a nitrogen overabondance by a 
factor of 7.44+4-08 in 2002 and 5 49 +4.51 in 20 03. We detected a 
modulation of the X-ray flux phased with the orbital motion, but 
could not find short term variability (from a few minutes to about 
one hour) that might have been expected from hydrodynamical 
instabilities of the wind-wind collision. 

In this paper, we discuss a set of high resolution optical spec- 
tra of Plaskett's star, in order to constrain the evolutionary sta- 
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Fig. 1. Normalized disentangled spectra of HD 47 129, over 3 distinct 
wavelength ranges. The secondary spectrum was vertically shifted by 
0.4 units for clarity. 



tus of this peculiar system, and to determine the role played 
by the Struve-Sahade effect. Indeed, a better knowledge of the 
processes that take place in HD 47129 would be useful to un- 
derstand the massive O+O binaries in general. The data have 
been obtained between 2002 and 2007 at the European Southern 
Observatory (ESO, Chile) and at the Observatoire de Haute- 
Provence (OHP, France). 

The paper is organized as follows. Section 2 describes the 
observations and data reduction. In Section 3, a radial veloc- 
ity curve, computed with a disentangling method, and the cor- 
responding orbital solution are presented. We discuss the spec- 
tral type of HD 47 129 in Section 4, the behavior of the emission 
lines using a Doppler tomography technique in Section 5 and the 
behavior of the absorption lines in Section 6. We use a model 
atmosphere code in Section 7 in order to fit theoretical spectra 
to the disentangled spectra of the components of Plaskett's star 
and discuss the evolutionary status of the system in Section 8. 
Finally, Section 9 summarizes the results and gives the conclu- 
sions. 



2. Observations and data reduction 

The journal of observations is presented in Table 1 . A part of the 
data was obtained by our team at ESO (La Silla) with the Fiber- 
fed Extended Range Optical Spectrograph (FEROS), mounted 
on the 1.5 m telescope (7 spectra in 2002, run ESO 068. D- 
095(A)) and 2.2 m telescope (5 spectra in 2006, run ESO 076.D- 
0294(A)). We have also added three FEROS archive spectra 
taken in 2006 (PI: Casassus, run ESO 076.C-043HA) and PI: Lo 
Curto, run ESO 076. C-0 164(A)). The wavelength domain of this 
echelle-spectrograph is [3550-9200] A. Typical exposure times 
range from 5 to 15 minutes, with a mean signal to noise ratio 
(SNR) equal to 150 in 2002 and 340 in 2006. The spectral re- 
solving power of FEROS is 48000. The data were reduced us- 
ing an improved version of the FEROS pipeline (see Sana et al. 
2003) working in the MIDAS environment. 
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The second part of the data was obtained in March 2007 at 
the OHP with the new echelle-spectrometer Sophie, mounted at 
the 1.93m telescope. The Sophie instrument features a 52.65 
grooves/mm R2 echelle grating and the spectra are projected 
onto an E2V 44-82 CCD detector. Thirty-nine spectral orders 
are currently extracted over a wavelength domain equal to 
[3872-6943] A, with a resolving power of 40000 in the High- 
Efficiency (HE) mode (Bouchy & the Sophie Team 2006). The 
mean SNR in these spectra is equal to 290. The data reduc- 
tion is done via an entirely automatic pipeline, adapted from the 
HARPS software designed at the Geneva Observatory. 

Finally, data were also obtained with the Aurelie spectrom- 
eter, mounted at the 1.52 m telescope at OHP and equipped 
with a 2048 x 1024 CCD EEV 42-20#3, with a pixel size of 
13.5 fim squared. All spectra were taken with a 600 1/mm grat- 
ing with a reciprocal dispersion of 16 A mm -1 , allowing us to 
achieve a spectral resolving power of 8000 in the blue range (6 
spectra in the [4450-4900] A domain), 9500 in the green range 
(7 spectra in the [5480-5930] A domain) and 11000 in the red 
range (1 spectrum in the [6340-6780] A domain). The spectra 
have a mean SNR equal to 300. The data reduction has been done 
with the MIDAS software as described in Rauw & De Becker 
(2004). Because of their lower resolution, the Aurelie data were 
not used for the disentangling and Doppler tomography tech- 
niques applied in this paper (see Sections 3 and 5). However, 
they provide a set of spectra well distributed in phase and we 
used them in order to check the accuracy of the results found 
with the FEROS and Sophie spectra. 



3. Orbital solution 

As for the main sequence early-type binary systems studied in 
Linder et al. (2007), the radial velocities (RVs) were determined 
in an iterative way using a disentangling algorithm (Gonzalez & 
Levato 2006) applied to three different wavelength domains of 
the FEROS and Sophie spectra ([4000-4220], [4360-4600] and 
[4900-5060] A, presented in Fig. 1 in which the most impor- 
tant lines are labeled). We choose these domains because they 
present a lot of helium lines that are also useful for the spectral 
type determination, and exclude complicated emission lines such 
as He ii A 4686 and Ha (see Sect. 5). For each of these domains, 
we built two specific masks 1 for the cross-correlation with the 
primary and secondary spectrum, respectively. This mask is dif- 
ferent for the primary and the secondary component because 
the nitrogen lines that are prominent in the primary spectrum 
are quasi absent from the secondary spectrum. The lines used 
for the RV calculation are Hd; Hei AA 4388, 4471, 4921, 5016; 
Hen AA 4200, 4542 and Si iv A 4089. For the primary star, we 
also added Hei AA 4009, 4026, 4120, 5048; Siiv A 4116 and 
Nm A 4097, 4196. The disentangling code converged to final 
separated spectra for each component and for each wavelength 
domain, and to corresponding radial velocities. 

The primary RVs show very little scatter between the three 
domains and the dispersion on individual data points are on aver- 
age less than 2.0 km s _1 . For the secondary, the situation is unfor- 
tunately less clear and the dispersion is significantly larger (most 
of the time of the order of 20km s~'). For three data points, the 
error exceeded 100 km s -1 and they have been completely ex- 
cluded from the analysis. There are several reasons for this sit- 
uation. One is the existence of changes in the line intensities 
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Fig. 2. Radial velocity curve of HD 47 129. Primary data are represented 
by filled symbols, while secondary data are represented by open sym- 
bols. Circles (resp. squares) stand for FEROS data (resp. Sophie data). 
The RVs have been determined through a disentangling of the primary 
and secondary spectra using the method of Gonzalez & Levato (2006). 
O = corresponds to the phase of the conjunction with the primary 
being in front. 

Table 2. Up: Orbital period of HD 47 129 calculated with the Fourier 
technique with all data available in the literature. Middle: Orbital so- 
lution of HD 47 129 computed from the RVs of the primary, obtained 
from the disentangling method. To refers to the time of the primary con- 
junction (O = 0). y, K and a sin i are respectively the systemic velocity, 
the semi-amplitude of the radial velocity curve, and the projected sepa- 
ration between the centre of the star and the centre of mass of the binary 
system. Quoted uncertainties are the l<x error bars. Below: Derived el- 
ements of the secondary orbit (see text). 



Parameter 


P (days) 


14.396257 


± 


0.000953 


e 







(fixed) 


To (HJD) 


2452331.547 


± 


0.028 


7p(kms _1 ) 


30.6 


± 


1.8 


K v (km s~') 


202.2 


+ 


2.6 


a P sin i (R ) 


57.5 


+ 


0.7 


F mass , P (Mq) 


12.3 


± 


0.5 


rms (km s" 1 ) 




8.6 




ys (kms -1 ) 


-6.0 


± 


9.6 


K s (km s- 1 ) 


192.4 


± 


6.7 


M s sin 3 i (km s" 1 ) 


47.3 


± 


0.3 


M P 


1.05 


± 


0.05 



1 These masks are synthetic spectra equal to except at the theoreti- 
cal wavelengths of the lines where they are equal to 1. 



of both the primary and secondary spectral lines. In particular, 
this leads to an imperfect subtraction of the spectral signature 
of the primary from the cross-correlation function, thus shift- 
ing the measured RV of the secondary toward the primary's RV. 
To circumvent this problem, we have scaled the primary spec- 
trum by the ratio between the height of the main peak in the 
cross-correlation function evaluated at individual phases and the 
equivalent in the cross-correlation function of the mean primary 
spectrum before subtracting the latter from the observed spec- 
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trum at the individual phases. This solves the issue for the pri- 
mary line variability, but the line variations of the secondary still 
affect the resulting RVs (Fig. 2). Another reason for larger dis- 
persion of RV's in the secondary star is the signifigant rotational 
broadening of the secondarys absorption lines, which makes it 
difficult to determine the centroi'd of the line. Table 1 shows the 
averaved RV values from the three domains for both the primary 
and the secondary star. These values were then used one more 
time as a fixed input to our disentangling code to derive a final 
set of self-consistent disentangled spectra presented in Fig. 1. 

The second step consisted in a refining the orbital period by 
considering all the data available in the literature (Struve 1948; 
Abhyankar 1959; Hutchings & Cowley 1976; Stickland 1987; 
Wiggs & Gies 1992; Bagnuolo & Barry 1996) in addition to the 
spectra presented in this paper. The new orbital period, equal to 
14.396257 ± 0.000953 days, has been determined with the gen- 
eralized Fourier technique (Heck et al. 1985; Gosset et al. 2001), 
especially adapted to unevenly spaced data. The error quoted 
here is empirical, and is equal to one tenth of the width of the 
peaks of the power spectrum. The theoretical formula of Lucy 
& Sweeney (1971) gives a smaller error value based on the sig- 
nal to noise ratio but actually, with the inhomogeneous data set 
considered here, the empirical error is preferred. 

The orbital solution was calculated with the program LOSP 
(Liege Orbital Solution Package, see Rauw et al. 2000), keep- 
ing the orbital period fixed at the above value. This program is 
a modified version of the algorithm developed by Wolfe et al. 
(1967). The calculation was made on primary data only, as for 
an SB1 system and the solution, given in Table 2, is close to the 
results found by Stickland (1987) with all the available data, and 
by Stickland (1997) with only the 1UE data. A major difference 
between their solution and ours is the systemic velocity y, which 
is ~ 8 km s _1 larger in our orbital solution. Our mass function is 
also higher than the mass function calculated with the IUE data. 

From the SB 1 solution, we have fitted the secondary radial 
velocities to the following relation: V rj s(0 = CiV r ,p(0 + c 2 where 
V r> P and V r? s are respectively the primary and secondary radial 
velocity as a function of time, c\ = — -^- with Mp, Ms the masses 
of the primary and secondary component respectively. We then 
found Q = H = 1.05 ± 0.05, ^ s = 192.4 ± 6.7 km s" 1 and y s = 

-6.0 ± 9.6 km s _1 . These results give M P sin 3 i - 45.4 ± 2.4 M 
and Ms sin 3 i = 47.3 ± 0.3 M . If we suppose an orbital inclina- 
tion between 69.3 and 72.7°(Bagnuolo et al. 1992), these num- 
bers lead to masses between 52 and 55 M Q for the primary com- 
ponent, and between 54 and 58 M Q for the secondary compo- 
nent. These results are summarized in the second part of Table 2. 

Orbital solutions of individual spectral lines, based on di- 
rect measures performed on the observed spectra, have also been 
investigated and no significant differences have been noticed. 
These measures were only possible on the primary component 
of the system, the secondary lines being too faint to obtain reli- 
able results with this method. 

The projected rotational velocities (v sin i) were derived 
using the Fourier technique (Simon-Diaz & Herrero 2007; 
Gray 2005). For the primary, we found that the width of the 
Hei AA 4471, 4921 and 5016 lines indicates v sin i = 75, 62 and 
60 km s _1 , respectively. The He n AA 4200, 4542 lines cannot be 
used here because of the severe blending with the nearby N in 
lines. As far as the secondary star is concerned, the Hei A 4471, 
4921 and Hen AA 4200, 4542 lines yield v sin/ = 245, 230, 
310 and 305 km s -1 respectively. Here, the latter values (Hen) 
are likely to be the correct ones because the spectrum of the 
secondary displays some emission in the wings of the Hei 



Table 3. Projected rotational velocities, derived using the Fourier tech- 
nique for some He I and He II lines. The errors have been estimated 
from the uncertainties on the position of the first zero in the Fourier 
transform of the observed profile. 





v sin 


(km s ' ) 




Primary 


Secondary 


He I ^4471 


75 ±5 


245 ± 15 


He I ^4921 


62 ±3 


230 ± 15 


He 1 ,15016 


60 ±3 


- 


He II A 4200 


- 


310 ±20 


He II A 4542 


- 


305 ± 20 



lines which leads to apparently narrower He i absorptions (see 
Table 3). 

4. Spectral type determination 

A determination of the HD 47 129 spectral type can be made 
following Conti's criterion (Conti & Alschuler 1971, and subse- 
quent articles), which is based on the ratio between the equiv- 
alent width (EW) of Hei A 4471 and Hen A 4542. We first 
measured the EWs directly on the observed spectra, after a 
standard deblending of both component lines by fitting two 
Gaussian profiles to the data when the stars are near quadra- 
ture phases. This leads to an 09 spectral type for the primary 
star and 06.5 for the secondary star. The uncertainties on these 
numbers correspond to a range for the primary between 08.5 
and 09. On the other hand, the possible interval of spectral 
types for the secondary star is larger, between 06 and 07.5. 
This is mainly due to the presence of nitrogen lines which 
complicate the secondary component measure of He n A 4542. 
However, it is relatively close to the secondary's spectral type 
determined by Bagnuolo et al. (1992) with their tomographic 
separation method, which is not the case for the primary compo- 
nent (they found 07.5 I + 06 I). The behavior of the He n A 4686 
line suggests that at least one component of the system is not 
a main sequence star (Fig. 3). The narrow emission line in 
the Hen A 4686 complex seems indeed to move along with 
the primary star and would thus indicate a luminous giant or 
supergiant. We applied Conti's 07-09.5 luminosity criterion 
(log W" = log EW(Siiv A 4089) - log EW(Hei A 4143)) to 
the primary measures of EWs and it indeed yields a giant (III) 
luminosity class. 

Furthermore, we applied the same criteria on the disentan- 
gled spectra, and noticed substantial differences. Measures gave 
an 08 III spectral type for the primary spectrum and 07.5 I for 
the secondary spectrum. This is not surprising since Fig. 1 shows 
two very similar spectra, even though the secondary lines are 
heavily broadened. 

Comparison with the criteria from the atlas of Walborn 
& Fitzpatrick (1990) (Hei A 4771/Hen A 4542 and 
Hei A 4026/Hen A 4200 EW ratios) further confirms the 
08 spectral type of the primary star with a huge enhancement 
of the nitrogen lines. The ratio Hei A 4338/Hen A 4542 in the 
primary spectrum argues against a luminosity class brighter 
than III, however it is interesting to note that the emissions 
seen in S iv AA 4486, 4504 seem to be more typical of a 
supergiant 2 . Most criteria thus indicate 08 III/I for the primary 
star. Concerning the secondary star, the He i A All 1/He n A 4542 



2 We note however, that the "selective" nature of these lines is not 
fully established (Walborn 2001; Werner & Rauch 2001). 
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ratio in the disentangled spectra clearly shows a spectral 
type later than 07, while the absence of S iv emission lines 
suggests that it has a luminosity class III. Finally, the secondary 
should be of 07.5 spectral type, rather than 08, because of the 
absence of Hei AA 4143, 4388. In conclusion, we favor the 
following spectral types: 08 III/I for the primary, 07.5 III for 
the secondary. 

We derived an optical brightness ratio, by comparing the 
EW ratio between the primary and secondary disentangled 
spectra and the same ratio evaluated for a sample of presum- 
ably single stars of the same spectral type (Conti & Alschuler 
1971; Conti 1973). The comparison was performed for the 
He i AA 4026, 4388, 4471 and Hen AA 4200, 4542 lines, and 
we found -^ = 1.9 ± 0.1. The primary star is thus brighter than 
the secondary star in the optical domain. This result is similar to 
the optical brightness ratio found by Bagnuolo & Barry (1996) 
with their tomography algorithm. 

From the average of the UBV photometric data compiled by 
Reed (2005), we obtain V = 6.06 ±0.01 andfi-V = 0.04 ±0.03. 
Assuming (B - V) = -0.27 (Martins & Plez 2006), we in- 
fer E(B - V) = 0.31 ± 0.03, which implies A v = 0.96 ± 0.09 
if we adopt Ry - 3.1. Assuming HD 47 129 to be a member 
of the Mon OB2 association, the distance modulus should be 
about 10.9 (Humphreys 1978). Hence, the absolute magnitude 
of HD 47 129 would be M v = -5.80 ± 0.09. From the optical 
brightness ratio inferred above, we then arrive at absolute mag- 
nitudes of M v ,p = -5.342 ± 0.105 and M v ,s = -4.645 ± 0.131 
for the primary and secondary respectively. The value for the pri- 
mary is somewhat (0.13 magnitude) fainter than typical 08 III 
V-band luminosities (Martins & Plez 2006) whilst the M v of 
the secondary is about 0.9 magnitude fainter than the typical 
07.5 III magnitude and would be much more typical of a main 
sequence or subgiant luminosity class. A better agreement be- 
tween the absolute magnitude and the spectral type can be ob- 
tained, at least for the primary star, by increasing the distance 
of the system. Indeed, if we suppose that the primary star has 
an absolute magnitude equal to a typical value of an 08 III star 
(Af v ,p = -5.47, Martins & Plez 2006), the distance modulus 
would be equal to 1 1 .0, which is still compatible with the dis- 
tance determined for other stars belonging to the Mon OB 2 
association (i.e. HD 46 202 and HD 47240, Humphreys 1978). 
Taking this new distance into account, the absolute magnitude 
of the secondary star is now equal to -4.72, which is still too 
faint for an 07.5 III star and is more typical of an 06.5-07 V 
star. If we adopt bolometric corrections of -3.16 and -3.25 for 
the primary and secondary respectively (Martins & Plez 2006), 
we infer bolometric luminosities of log L/L G = 5.35 ± 0.04 and 
5.09 ± 0.04 respectively. 

At this stage, it is worth emphasizing that the rapid rotation 
of the secondary (see Section 3) could lead to a rotational flatten- 
ing that would make this star appear hotter near the poles than 
at the equator. As a result, the spectral type determined from 
the He i/He n line intensity ratios might actually be biased. We 
applied the equations given by Maeder (1999, Appendix A), as- 
suming that the inclination of the rotation axis is the same as the 
orbital inclination (i.e. 71°), and that the angular velocity is con- 
stant. We then found that, for an equator radius equal to 14.1 R 
(which is a typical value of an 07.5 III star, Martins et al. 2005), 
the polar radius is equal to 13.2 R Q . The von Zeipel theorem ap- 
plied to our case indicates a difference in effective temperature 
equal to ~ 2 500 K between the pole and the equator (von Zeipel 
1924). This means that the temperature of the secondary's pole 
would be equal to 35 500 K (since we derive T e ff = 33 000 K, 
see Sect. 7), which is compatible with an 07 V-07.5 V star. 



The secondary star of HD 47 129 could thus be an 07 V star 
that shows an apparent 07.5 III spectrum because of its rapid 
rotation. 

In conclusion, the discrepancy between the spectral type of 
the primary star and its absolute luminosity can be solved by 
adopting a distance modulus of 1 1.0 for the system. For the sec- 
ondary star this solution is not sufficient, but the rest of the dis- 
crepancy could be due to the large rotational velocity of the star, 
that would induce a strong gravity darkening leading to a lower 
temperature (and hence apparently later spectral type) over the 
equatorial regions. The same effect also yields a lower gravity 
near the equator than close to the poles, hence leading to an ap- 
parent giant luminosity class. 

Nevertheless, the rather modest luminosities are still surpris- 
ing given the minimum dynamical masses of the stars inferred 
from our orbital solution (45 and 47 M Q , see Sect. 3): both stars 
would appear underluminous for their dynamical masses. The 
best estimate of the inclination is 70°(Bagnuolo et al. 1992), 
leading to masses of 55 and 57 M Q , and thus enhancing the 
discrepancy. Actually, this problem can be formulated in a dif- 
ferent way: the stars have spectral types that are too late for 
their dynamical masses, which are more typical of giant or main 
sequence 03 stars. Increasing the distance would imply that 
Plaskett's star does not belong to the Mon OB2 association, 
which is unlikely. Another still unknown effect might thus be 
present here. 

5. The emission lines and Doppler tomography 

HD 47 129 displays prominent phase-locked profile variability 
of the strongest emission lines (Ha, Hen A 4686, Hei A 5876) 
in its optical spectrum. As becomes clear from a montage of 
the observed profiles (Figs. 3 and 4), the Ha and He n A 4686 
lines consist of several components: a narrow emission line that 
roughly follows the orbital motion of the primary star, and two 
broad emissions (a blue-shifted and a red-shifted one) that do 
not seem to move with the orbital cycle. 

Further insight into the origin of these variations can be ob- 
tained from the Doppler tomography technique which allows 
to map the formation region of an emission line in velocity 
space. Here we use an implementation of the Doppler tomogra- 
phy based on a Fourier filtered back projection algorithm (Home 
1991) already used by Rauw et al. (2002, 2005) for the analysis 
of line profile variations of HDE 228 766 and WR 20a. For this 
purpose, we adopt a reference frame centered on the centre of 
mass of the binary with the x-axis pointing from the primary to 
the secondary and the positive y-axis pointing along the direc- 
tion of the secondary's orbital motion. The Doppler tomography 
technique assumes that the phase dependence of the radial veloc- 
ity v(<t>) of any gas flow that is stationary in the rotating frame of 
reference of the binary can be described by a so-called 'S-wave' 
relation: 



v((J)) = Vx cos (2 n <t>) - v y sin (2 n O) + v z 



(1) 



where <t> is the orbital phase (<J> = at conjunction with the 
primary being in front), whilst v x and v v reflect the projected ve- 
locity components along the x and y axes and v z corresponds to 
the apparent systemic velocity of the line under investigation. A 
Doppler map consists of a projection of the S-wave relation of 
Eq. 1 onto the (v x , v y ) plane (for which v z = y). Each pixel in 
a Doppler map, specified by its velocity coordinates is associ- 
ated with a particular S-wave (see e.g. Home 1991 for a detailed 
discussion of the method). 
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Fig. 3. Hen A 4686 line in the spectra of HD 47 129 from FEROS (F) 
and Sophie (S) observations. The orbital phases are indicated to the 
right of each panel. 
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Fig. 4. Same as Fig 3, but for the Ha line in the spectra of HD 47 129. 

To build the Doppler maps of the Hen A 4686 and Ha 
emission lines (Fig. 5), the data were weighted so as to ac- 
count for the uneven sampling of the orbital cycle. The Doppler 
maps of both lines present their highest peak near the centre 
of mass of the primary with a slightly negative v x (-60.0 and 
-68.0 km s 1 for Hen A 4686 and Ha respectively) and v y very 
close to the amplitude Kp of the primary's radial velocity curve 
(202.5 km s _1 ). This peak is associated with the narrow emis- 
sion component visible in both lines and discussed above. In 
addition, both Doppler maps reveal a roughly annular region 



with patchy emission structures. This annulus is most promi- 
nent for the Hen A 4686 line. The structure is centred on a 
point with coordinates (v x , v y ) = (0, +25) km s -1 , with inner and 
outer radii of about 400 and 600 km s respectively. This fea- 
ture is responsible for the two broad, roughly stationary emis- 
sions discussed above. In the Ha Doppler map, the structure ap- 
pears more elongated and less well defined. In addition, the Ha 
Doppler map shows an extended 'emission bridge' connecting 
the peak at (v x , v v ) = (-68, -198) km s _1 with the annulus near 
(v x , v y ) ~ (-450, +450) km s _1 . 

Let us now turn to the possible interpretations of the vari- 
ous features seen in the velocity maps. First of all, we consider 
the peaks of the maps that occur at slightly negative v x with 
respect to the motion of the primary star. The velocity coordi- 
nates of this emission region hence indicate a motion from the 
secondary toward the primary. A possible interpretation could 
thus be that this emission arises from a region where the pri- 
mary wind is compressed by the interaction with the secondary 
wind (which is flowing toward the primary). In terms of a col- 
liding wind interaction, this would correspond to the head of the 
shock which should indeed be displaced toward the primary. The 
narrow emission line component thus likely forms near the re- 
gion where most of the X-rays are emitted (Linder et al. 2006). 
This scenario is in agreement with the one proposed by Wiggs 
& Gies (1992). Measures of radial velocities performed on the 
FEROS, Sophie and Aurelie spectra also showed that the nar- 
row emission lines come from a region of the wind close to the 
primary star. Actually, both S-waves are similar to the primary's 
RV curve presented in Fig. 2, but with a lag in phase of 0.04 for 
He ii A 4686 and 0.06 for Ha. This result corresponds to the col- 
liding wind model of Wiggs & Gies (1992) but with an emission 
zone of Ha maybe a little closer to the system's centre of mass 
than the emission zone of He n A 4686. 

An alternative explanation for the narrow component could 
be a 'hot spot' due to the impact of a gas stream from a Roche 
lobe filling component on its companion. However, HD 47 129 
is currently likely in a post Roche lobe overflow (RLOF) config- 
uration where the primary transferred material to the secondary 
(Bagnuolo et al. 1992). The (v x ,v y ) components of the narrow 
emissions are therefore at odds with such a scenario. 

A striking feature is the annular-like structure seen in the 
Doppler map that results from the roughly stationary double- 
peaked feature in the spectrum. Actually, both peaks display 
some variation in intensity and radial velocity, but they cannot 
be directly linked to the orbital phase of the system (this phe- 
nomenon is also observed in the absorption lines, see Sect. 6). 
This kind of structure could reflect the existence of an accretion 
disk around one of the stars. Indeed, it is believed that Plaskett's 
star has undergone mass exchange through a RLOF of the pri- 
mary component in the past. However, it seems unlikely that this 
interaction would still be ongoing, especially since several fea- 
tures rather hint at material flowing from the secondary toward 
the primary (see above), i.e. in the opposite sense compared to 
what is expected for the RLOF of the primary. Another point 
concerns the fact that the annular structure does not appear to 
be centered on any of the two stars. A more likely explanation 
could be that the secondary star has a wind flattened by its ex- 
tremely high rotational velocity and hence compressed into its 
equatorial plane in a way similar to the 09.5 V star HD 93 521 
(Bjorkman et al. 1994). The outer regions of such a flattened 
wind, where the Ha line forms, would probably be heavily dis- 
turbed by the presence of the primary star orbiting around the 
secondary. However, the inner regions of the wind, where the 
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Fig. 5. Top: Doppler map of the Ho- emission line in the spectrum of 
HD 47 129. The crosses correspond to the radial velocity amplitudes of 
the centre of mass of the primary and secondary. The shape of the Roche 
lobe in velocity-space (thick dashed line) has been calculated for a mass 
ratio (secondary/primary) of 1.05. The Doppler map was computed with 
v z set to +30 km s~' . The black dot indicates the position of the highest 
peak and the contours correspond to levels of 0.95, 0.80, 0.65, 0.50, 
0.35 and 0.20 times the maximum emissivity. Bottom: Same, but for 
the He ii A 4686 emission. 




Fig. 6. Schematic pole-on view of the interaction in HD 47 129 (not to 
scale). The sense of the motion of the two stars is given by the arrows. 
The secondary star has a rotationally flattened wind. The He n A 4686 
and Ha emission regions are highlighted and the wake due to the mo- 
tion of the primary across the secondary wind is indicated by the thick 
curve. The X-ray emitting region near the surface of the primary is also 
indicated (undulated arrows). 



As we have seen in Section 4, the primary star is intrinsi- 
cally more luminous than the secondary star and should hence 
have the stronger wind. Furthermore, the analysis of ROSAT- 
HRI archive data and of XMM Newton data have shown that a 
decrease in the X-ray flux seems to occur when the primary star 
is in front, which also indicates a stronger primary's wind, ab- 
sorbing more the X-ray flux than the secondary's wind (Linder 
et al. 2006). On the other hand, both the estimation of the com- 
ponent's mass loss rate in Linder et al. (2006) and the optical 
emission lines suggest that the secondary star has the stronger 
wind and thus that the wind interaction zone is displaced toward 
the primary star. This apparent contradiction can be solved eas- 
ily if the bulk of the secondary's wind is indeed concentrated in 
the equatorial plane as suggested above. 

Finally, the bridge-like structure seen in the Ha map is most 
likely produced by the wake created by the primary star moving 
across the flattened wind of the secondary star. 

A schematic view of a possible wind interaction between the 
components of Plaskett's Star, where the secondary wind is flat- 
tened by the action of its rapid rotation, is shown in Fig. 6. 



He n A 4686 line forms, would be less affected, thus explaining 
why the feature is best seen in the Doppler map of the latter 
line. We performed numerical computations that showed that it 
should be possible to qualitatively reproduce this annulus with 
a combination of the wind's rotation velocity and the wind's ex- 
pansion velocity. However, it was not possible to obtain exactly 
the same annulus because the variation of the wind's rotation 
velocity with the distance from the secondary star as well as the 
influence of the UV primary radiation, which is supposed to re- 
duce the secondary wind velocity, were unknown. 



6. The absorption lines and the Struve-Sahade 
effect 

The S-S effect is known to appear in massive early type O + O bi- 
naries and impacts their fundamental parameters, such as masses 
and luminosity ratios. A systematic study of this particular effect 
has been made by our team on main sequence systems (Linder 
et al. 2007), and we have determined that the S-S effect did not 
affect all studied lines of the spectra and that in most cases both 
stars of the system were concerned by the phenomenon. In or- 
der to continue the study, we analyse here the optical spectra of 
Plaskett's star which is one of the giant/supergiant O-type bina- 
ries that have been reported to display the Struve-Sahade effect. 
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Fig. 7. Observed FEROS (a, b and c) and Sophie (d) helium lines of Plaskett's star (solid line) and disentangled spectra appropriately shifted in 
phase (red dotted lines). The lower panels yield the residuals (observed-reconstructed) in units of the normalized spectrum. 



Measures of equivalent widths made on all observed spec- 
tra (with FEROS, Aurelie and Sophie) showed first that, even 
with high resolution and reasonable S/N ratio data, separating 
the primary and secondary contribution to a given line is not 
a simple task. Actually, the secondary spectrum is so broad- 
ened and sometimes so faint that the measures show a large 
dispersion. The method applied in Linder et al. (2007) was 
thus not optimal to check for the presence of the S-S effect in 
HD 47 129. Indeed, the standard deblending of both compo- 
nents by fitting two gaussians on the data showed that, among 
all tested lines (Hei AA 4026, 4388, 4471, 4713, 4921, 5016; 
Hen AA 4200, 4542 and Nm AA 3999, 4003, 4195, 4379), only 
two absorption lines seemed to be slightly stronger at O = 
0.75 (when the secondary star is receding): Hei A 4026 and 
He i A 4388. 

On the other hand, the disentangled spectra can be consid- 
ered as a mean of the observed spectra. We can thus check for 
the presence of the S-S effect by comparing these mean pro- 
files to the observed ones for each specific orbital phase. For this 
purpose, we have reconstructed the spectra by shifting the disen- 
tangled profiles of the primary and secondary by the appropriate 
radial velocities of these components. In this way, we find that at 
specific orbital phases, the observed line profiles often differ sig- 
nificantly from the reconstructed ones. The strongest variations 
are clearly associated with the secondary star's lines that vary 
both in width and intensity. However, to a lesser extent, some 
variations are also observed for the primary star's lines (see e.g. 
the N in A 4197 line in Fig. 7). It is not completely clear to what 
degree the variations of the absorption lines are phase-locked. 
Indeed, whilst some observations taken at close phases exhibit 
similar deviations from the reconstructed profiles, other obser- 
vations with almost identical phases reveal totally different pro- 
files. As a result, the line profile variations appear dominated by 
a modulation that is either stochastic or at least not tied to the 
orbital phase. 

As stated above, the variations mostly affect the lines of the 
secondary star and the strongest deviations are seen in the wings 
of the Hei lines. On some spectra, the observed profile appears 
broader than the average, whilst on other data the reverse situa- 
tion is observed. There are a number of possible scenarios to ex- 
plain variations that are not linked to the orbital cycle. The fact 
that the variability mostly affects the wings of the secondary's 
He i lines could be a hint that it actually stems from variations 
of weak emission wings. Indeed, such wings are clearly seen to 
flank some of the stronger Hei lines (see Fig. 1). In the context 
of our suggestion that the secondary of HD 47 129 is surrounded 
by a flattened wind, the variability of these emission wings could 



then reveal a rather complex velocity structure of this flattened 
outflow. If correct, this interpretation would imply that the appar- 
ent variations of the absorption lines in HD 47 129 actually result 
from the effect of the turbulence and other velocity perturbations 
on the emissions formed in the flattened wind of the secondary 
star. Alternatively, the variability could be tied to the rotational 
cycle of the secondary or to non-radial pulsations of the latter. 
Concerning the first possibility, we note that the uncertainties on 
the rotational period of the secondary prevent us from further ex- 
ploring this possibility. Non-radial pulsations in rapidly rotating 
O stars can induce prominent line profile variability (see e.g. the 
case of HD 93 521, e.g. Howarth & Reid 1993). However, the 
amplitude of the line profile variations in HD 47 129 are rather 
substantial (several percent of the combined continuum, com- 
pared to less than one percent for HD 93 521) and accounting 
for the brightness ratio between the components of Plaskett's star 
would imply intrinsic amplitudes even about three times larger. 



7. Modeling the spectra with an atmosphere code 

We used the model atmosphere code CMFGEN (Hillier & Miller 
1998) in order to derive the wind and photospheric properties of 
both components of HD 47 129 from the disentangled spectra. 
CMFGEN calculates non-LTE atmosphere models that include 
winds and the treatment of line-blanketing. Its main character- 
istics and hypothesis can be found in Martins et al. (2008). In 
addition to our optical spectra, we made a comparison between 
IUE archive spectra to which we have applied our disentangling 
method and theoretical spectra in order to determine the wind pa- 
rameters. We first performed the tests taking the luminosities we 
observed for the stars into account, and then luminosities equal 
to the average value for their spectral type (08 III/I + 07.5 III, 
Martins et al. 2005). The results provided by these two models 
are not significantly different, and thus fall within the same un- 
certainties. 

Table 4 gives the main results we obtained. The effective 
temperatures and the effective gravity are compatible with the 
values given by Martins et al. (2005) for the spectral types we 
found in Section 4. The IUE spectra allowed us to calculate the 
wind parameters of the primary star. The / parameter, or filling 
factor, is the volume occupied by matter in the wind divided by 
the total volume of the wind. Its value has been adopted in order 
to best reproduce the N iv A 1 720 and O v A 1371 lines, which are 
particulary sensitive to clumping. The mass loss rate we found 
here is three times weaker than the theoretical mass loss rates 
given by Linder et al. (2006). The uncertainty on the quantity 
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Table 4. Some parameters of Plaskett's star, determined from simula- 
tions with the atmosphere code CMFGEN. T e j stands for the effective 
temperature, M for the mass loss rate, / for the filling factor and v„ 
for the wind terminal velocity. The abundances of the listed elements 
are in number (relative to H), and compared to the solar abundances of 
Grevesse et al. (2007). 



Primary 



Secondary 



Spectral type 


08 III/I 


07.5 V/III 


log L/L e 


5.35 (fixed) 


5.09 (fixed) 


T eff (K) 


33500 ± 2000 


33000 ± 2000 


logg 


3.5 ±0.1 


3.5 ±0.1 


He/H 


0.10 ±0.03 


0.15 ±0.05 


He/He 


1.17 ±0.35 


1.76 ±0.53 


C/H 


5.0 x 10~ 5 ± 1.5 x 10~ 5 


2.45 x 10~ 4 (fixed) 


c/c Q 


0.2 ±0.1 


1.0 (fixed) 


N/H 


1.0 x 10- 3 ±0.3x 10~ 3 


1.0 x 10- 5 ±0.3x 10~ 5 


N/N 


16.6 ±5.0 


0.2 ±0.1 


Mg/H 


1.9 x 10~ 5 ±0.6x 10~ 5 


- 


Mg/M go 


0.6 ±0.2 


- 


M(M Q yr') 


1.7 x 10- 7 ±0.3x 10~ 7 


- 


/ 


0.05 ± 0.03 


- 


Voo (km s ') 


2450 


- 



log -M= is about 0.2-0.3 dex and was estimated as the change 

beyond which the corresponding synthetic spectrum does not fit 
anymore the observed spectrum. This can be considered as a 
3 sigma uncertainty. The abundances are listed by number rel- 
ative to hydrogen and the corresponding error bars are equal 
to ~ 30%. They have been determined by taking different val- 
ues of the abundance for each element and comparing the re- 
sults with the disentangled spectra. The nitrogen abundance we 
found in the primary star is equal to approximately 16 times 
the solar abundance (taking the solar abundances of Grevesse 
et al. 2007). The nitrogen lines used to obtain this result are 
Nin AA 4511, 4515, 4518, 4523, 4528. This nitrogen enhance- 
ment is consistent with the overabundance of about 6 obtained 
with XMM-Newton data for the global X-ray emitting plasma 
of Plaskett's star (Linder et al. 2006). Indeed, in the X-ray data 
it was not possible to separate the contribution of the primary 
and secondary plasma in the emission, and thus the overabun- 
dance reported by Linder et al. (2006) was for the entire system. 
The lines used to determine the carbon abundance in the primary 
star are Cm AA 4070, 4152, 4156, 4163. It appears that there is 
a clear carbon underabundance (3 % of the solar abundance) in 
this component. The N enrichment and C depletion of the pri- 
mary component are expected from CNO processing in massive 
stars (Meynet & Maeder 2003). In the secondary star, it is the 
helium abundance that shows an overabundance, which is sur- 
prisingly coupled with a large underabundance in nitrogen. The 
abundance in carbon has been fixed to the solar value, because 
there is no relevant line in the secondary spectrum that could 
serve for a determination. 

Figs. 8 and 9 show the comparison between the observed 
and the theoretical spectra in the optical and UV domain, respec- 
tively. The correspondance between the optical part of the the- 
oretical spectrum and the disentangled spectra is good for both 
components. For the IUE spectra, the agreement with the pri- 
mary component is quite good, but the secondary star spectrum 
is not correctly represented and the secondary's wind parame- 
ters could not be determined, due to the puzzling shape of the 
UV spectrum of the secondary, where the classical wind lines 
(N v A 1240, Si iv AA 1393, 1403, C iv A 1548, 155 1) are not seen. 
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Fig. 8. Comparison between the disentangled optical spectra of 
HD 47 129 (solid line) and the theoretical spectra (red dotted line) com- 
puted with the CMFGEN atmosphere code. The secondary spectrum 
was vertically shifted by 0.3 units for clarity. 
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Fig. 9. Same as Fig. 8 but for the UV domain. The disentangled spectra 
have been built from IUE archive spectra of HD 47 129. The secondary 
spectrum was vertically shifted by 1.5 units for clarity. 



The X-ray luminosity has been taken into account in the models, 
which allows a better fit of the N v A 1240 line (log — £■ = -7, 
Chlebowski & Garmany 1991) although the blue part of the line 
is contaminated by the Lyman a interstellar absorption line. 

8. The evolutionary status of HD 47 129 

Some of the most important properties of Plaskett's Star suggest 
that it is currently in a post Case A RLOF evolutionary stage. 
The first one is the high rotational velocity of the secondary star. 
In fact, when the component that fills up its Roche lobe trans- 
fers material to its companion, there is also a huge amount of 
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angular momentum that is simultaneously transferred. As a re- 
sult, the accretor is spun up, whilst the mass donor spins down 
(e.g. Wellstein 2001; Petrovic et al. 2005). In a wide binary sys- 
tem, the accretion efficiency can be rather low: only a fraction 
of the mass lost by the donor is actually accreted, but this is 
nevertheless sufficient to spin up the accretor to substantial ro- 
tational velocities. The accretion efficiency is also reduced by 
the fast rotation of the accretor. Indeed, accretion would fully 
stop if the accretor is spun up to its critical velocity (Petrovic 
et al. 2005). In conclusion, the secondary of HD 47 129 was 
likely spun up to its present-day high rotational velocity dur- 
ing a previous mass transfer event. The fact that v sin i remains 
high is an indication that either the mass transfer stopped only 
very recently, or the synchronisation of the rotation due to tidal 
interaction is very inefficient in this case. If the average accre- 
tion efficiency was indeed as low as 10%, then the present-day 
properties of HD 47 129 could only be explained if the primary 
was initially very massive (around 80 M Q ) and the secondary had 
about 60 M Q . However, it must be stressed that a system with 
as high a primary mass would be expected to undergo an LBV 
phase rather than a RLOF (Vanbeveren 1991). 

Second, the He overabundance of the secondary star is likely 
related to a previous RLOF event. Indeed, the secondary ap- 
pears less evolved than the primary in terms of CNO process- 
ing: while the latter clearly shows a N overabundance and a 
C depletion typical of an advanced state of chemical evolution 
(e.g. Meynet & Maeder 2003), the former appears rather un- 
evolved, at least in terms of N enrichment. On the other hand, 
the secondary component has a higher He content compared to 
the primary for which He/H has not yet changed from its ini- 
tial value. From this, one can conclude that the He enrichment 
of the secondary star is not the result of classical evolution, but 
must originate from a mass transfer. The only caveat in this sce- 
nario is the underabundance of nitrogen in the secondary which, 
at present, cannot be explained by a RLOF event. Hunter et al. 
(2008) discuss the surface abundances of a sample of early O- 
type stars in the LMC as a function of their projected rotational 
velocities. Whilst most of the objects showed a behavior consis- 
tent with the expected chemical enrichment of the surface lay- 
ers due to rotational mixing, Hunter et al. (2008) identify two 
groups of stars with peculiar properties. Their group 1 consists 
of O-type stars that are fast rotators but do not show a significant 
nitrogen enhancement in their spectra. On the contrary, group 2 
of Hunter et al. (2008) contains O-type stars that display a rather 
low v sin;, whilst simultaneously showing a significant nitro- 
gen enrichment. The properties of these objects appear some- 
what related to those of the components of Plaskett's star. The 
rotational velocity and the surface composition of the primary 
star of HD 47 129 are reminiscent of those of the group 2 stars, 
whereas the secondary component of HD 47 129 has properties 
that would place it among group 1 objects. Furthermore, Langer 
et al. (2008) noted that mass transfer can produce rapidly rotat- 
ing stars that are not strongly nitrogen-enriched if the transfer is 
highly non-conservative. 

9. Conclusion 

We have presented a spectroscopic analysis of the massive O- 
type binary HD 47 129. 

The radial velocities of both components have been inferred 
using a disentangling method. The orbital solution gave two 
massive components of 54 M and 56 M Q for the primary and 
the secondary star respectively (assuming an inclination of 71°). 
The projected rotational velocities, determined using a Fourier 



technique, are very different from each other (~ 75 km s -1 for 
the primary, and ~ 300 km s for the secondary). 

The spectral types, calculated from the two disentangled 
spectra, gave 08 III/I + 07.5 III, which are too late spectral types 
for the dynamical masses. This discrepancy cannot be solved by 
taking into account the large rotational velocity of the secondary 
star. Whilst this effect indeed acts on the secondary luminosity, 
it is not sufficient to explain such a difference and it explains 
nothing concerning the primary star. 

The optical spectra display a phase-locked variability in the 
strong emission lines He n A 4686 and Ha. A Doppler tomog- 
raphy technique allowed us to obtain the Doppler maps of both 
lines which present their highest peak near the primary compo- 
nent velocity. Its velocity coordinates indicate a motion from the 
secondary toward the primary, which can be interpreted by a col- 
liding wind interaction, where the shock is displaced toward the 
primary star. The Doppler map also shows an annular feature, 
almost centered on the system's center of mass, that seems to 
be another consequence of the secondary's rotation. Indeed, it 
likely represents the secondary's wind compressed into its equa- 
torial plane by the fast rotation. 

Furthermore, we used the high resolution spectra in order to 
study the Struve-Sahade effect in this system. Actually, while 
Plaskett's star presents strong variations in its absorption lines, 
it is difficult to connect them to the 'standard' definition of the 
Struve-Sahade effect, which is: "the apparent strengthening of 
the secondary spectrum of the binary when it is approaching, 
and the corresponding weakening of its lines when it is reced- 
ing" (Bagnuolo et al. 1999). The variations are visible in both 
stars, even though they are stronger in the secondary's spectrum, 
but it is not clear that they are phase-locked. They appear domi- 
nated by a modulation that is either stochastic or at least not tied 
to the orbital phase. If correct, the interpretation of a flattened 
secondary wind would imply that these apparent variations result 
from the effect of the turbulence and other velocity perturbations 
on the emissions formed in the disk around the secondary star. 
The deformation of the star itself due to the rotation can also be 
in this case an additional, but not determining, perturbation. In 
the case of the main sequence stars of Linder et al. (2007), the 
Struve-Sahade effect was already found to be due to an inhomo- 
geneous distribution of the line formation, because of the defor- 
mation of the star and/or mutual heating, in at least three out of 
four studied objects (HD 100213, HD 159 176 and DH Cep). 

Finally, we compared our observations with theoretical spec- 
tra computed with the model atmosphere code CMFGEN. We 
then determined effective temperature and chemical abundances 
for both components, and wind properties for the primary star. 
The nitrogen overabundance in the primary star, already detected 
in the XMM-Newton data, is confirmed in the optical spectra. 
This comes with an helium overabundance of the secondary star, 
which suggests a past RLOF in the system, from the primary star 
to its companion. 

In conclusion, we can say that most of the peculiarities in 
the spectra of HD 47 129 can be explained today by the evolu- 
tionary status of the system (a post-case A RLOF) and by the 
fact that the projected rotational velocity of the secondary star 
is quite important. However, the large dynamical masses of both 
components are still not compatible with their spectral types and 
luminosities. 
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